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A microscopic model for thermal excitation of vibrational ground 
state of a molecule interacting with a condensed medium is de- 
veloped. The master equation for evolution of occupancies of the 
vibrational levels is derived. The rate constant limiting the pro- 
cess of molecular thermal excitation is determined analytically. It 
is shown that while this quantity is independent of temperature at 
low temperatures coinciding with a rate speed limit for quantum 
transitions due to uncertainty principle, at high temperatures it 
rather linearly rises with temperature according to the Einstein's 
relation for adiabatic transitions. 



Introduction 

Considering kinetic and dynamic processes in differ- 
ent atomic and molecular nanostructures, for example, 
in molecular impurity centers, metal or semiconductor 
nanoparticles, biological macromolecules etc., requires 
the knowledge of microscopic mechanisms for vibrational 
relaxation of the system at limiting stages of its evolu- 
tion. In recent years this problem became very impor- 
tant, particularly in the context of invention and impli- 
cations of superfast heat generation in the close vicinity 
of gold nanoparticles under impulsive irradiation [1, 2]. 
From the physics point of view, such process is analo- 
gous to instantaneous forming of a point source of heat 
energy for molecules [3] followed by their thermal excita- 
tion and subsequent transition to higher vibrational lev- 
els. Therefore, in the frame of microscopic approach one 
has to solve the problem of impact of thermodynamic 
fluctuations in condensed medium which is at certain 
temperature on the process of vibrational relaxation of 
molecule, that was in ground vibrational state initially. 

Many attempts were made to developing microscopic 
models for kinetic description of vibrational relaxation 
in molecules [4-9]. However, attention was mainly paid 
to analysis of intra- and intermolecular anharmonic- 
ity [5, 6, 9, 10] and role of spectral factors [9-12] in- 
volved in fluctuation-dissipation relations [9-11, 13]. At 



the same time, some specific types of random motions, 
particularly thermal noise (that is nonlinear over the 
molecular degrees of freedom), do not obey fluctuation- 
dissipation relations in the general case [14]. To ap- 
proach such types of motion, in this communication 
we use a concept of generalized open quantum sys- 
tem that interacts with a heat bath [15, 16]. Within 
this framework one can take into account both bilin- 
ear intermolecular anharmonicity (that, according to the 
fluctuation-dissipation theorem, causes one-phonon re- 
laxation transitions between molecular vibrational lev- 
els) and thermal noise (that randomly shifts these lev- 
els). We construct the microscopic model for thermal 
excitation of molecular vibrational ground state and de- 
rive the master equation for level occupancies. More- 
over, we determine the rate constant of thermal exci- 
tation analytically and show that, in the general case, 
this quantity exhibits simple temperature dependence 
- {2a;o + {2u}c/7T)exp{huJc/2kT)[exp{huJc/kT) - l]'^}, 
where h and k denote the Planck and Boltzmann con- 
stants, respectively, while the frequencies of vibrations 
in given molecule ujq and molecular environment ujc are 
determined only by structural factors and do not depend 
on the temperature T. 

1. Hamiltonian of a system 

We deflne a ground vibrational state as the state 
with highest kinetic energy K{q) = — ^^^r and 
lowest potential energy. The latter in harmonic 
approximation takes the form U{q,{Q}) = Uq + 
(l/2){a2(7(<z,{Q})/9g2}l,=,o(9 - 9o)^ where Uo = 
Uiq,{Q})\q=qQ — 0. Here q and qo denote the nor- 
mal nuclear coordinate of a given molecule (for simplic- 
ity we consider only one such coordinate with reduced 
mass m) and its equilibrium value found from the con- 
dition {■§^U{q,{Q})}\q^qg — 0, respectively. Potential 
energy depends also on displacements of the normal nu- 
clear coordinates of molecules in the environment {Q} 
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from their equihbrium positions {Qo}, that is necessary 
for occurrence of relaxation transitions between molec- 
ular vibrational levels. The interaction part of expan- 
sion of potential energy over mutual displacements (q — 
'?o)({Q} — {Qo}) responsible for such transitions appears 
in the form = { g J{q} U (g, {Q})}|g=qo.{Q}={Qo} (g ' 
Q'o)({Q} — {Qo}) [5, 6]. In this expansion the term 
~ (g— go)^ is excluded, that may be valid when the center 
of molecular vibrations coincides with the center of in- 
version. Also in Vrei there are absent terms, that do not 
make a key contribution to kinetics of relaxation transi- 
tions on the limiting stage of establishment of the final 
equilibrium, as well as terms, that can be taken into ac- 
count by renormalization of the Hamiltonian of the envi- 
ronment ii{{Q\) [6] (i.e., terms (g-<?o)^({(3}-{Qo}); 
(g-Qo)* and ({Q}-{go}); (<7-go)({Q}-{go})' respec- 
tively). Thus, in the representation of oscillator wave 
functions \n) (where n = 0,1,2,... is the number of the 
vibrational state) the respective Hamiltonian of entire 
system (molecule+environment+interaction) takes the 
form: 

H{t)=H^{t)+HT + Vrel (1) 

In this expression 

Ho{t)^^hL^{t){n + l/2)\n){n\ (2) 

n 

represents the main Hamiltonian of a molecule under 
consideration which is a stochastic operator depending 
on the random realizations of molecular movement in 
the environment along quasiclassical nuclear coordinates 
{Q{t)} involved in the expression for frequency of molec- 
ular vibrations 

[Lo{t)Y ^ [io{Q{tW - ^{^C/(g,Q(0)}l9=9o (3) 
Besides, in equation ([T]) 

HT = Y,hnx{h+bx + 1/2) (4) 

A 

is the Hamiltonian of a condensed medium, that in the 
harmonic approximation can be represented as a heat 
bath having an infinite set of the normal vibrations 
(phonons) with frequencies Quantities &J and h\ 
refer to the operators of creation and annihilation of the 
respective Ath mode. The last term in JT]) 

K-e; =^%AV^^(|n)(n+l| + |n+l)(n|)(6++6A) (5) 

nX 



represents the operator of relaxation transitions in the 
one-phonon approximation, that describes transitions 
between the nearest vibrational levels. The correspond- 
ing parameter of bilinear intermolecular anharmonicity 



X\ = V^^^o^^A (6) 

characterizes the dependence of relaxation transitions 
on both the frequencies of normal vibrations in a 
heat bath and the shifted frequency of molecular 
vibrations wq = a;((5(t))|{Q}={Q(,}. This frequency shift 
appears due to interaction of a molecule with a random 
mean field of thermal motions in the environment. The 
respective non-shifted frequency is wo = vjq — a, where 

^0 = T^{w^^°Hq)}U=<io with [/(0)(g) being the time- 
independent part of the potential energy. The latter can 
be found from the condition { ^^^^ U {q, {Q})}\q=qo = 0- 
The corresponding standard mean field deviation 

is r^'^ = ;fe{\/[^^(9'{^?W})]l«=^o,{Q}={Qo} " 

^[^C/(o)(g)] while 6 

1 ({3"C/(g,{Q})/9gaQ;,}|,^,„.{Q}^{Q„})^ 

4 {d''U{q,{Q})/dq''}\,=,oid'Hi{Q})/dQl}\iQy^^Q^y ^ ^ 

fines the so called parameter of nonadiabaticity. The 
latter couples the molecular vibrations with those of 
the environment (in the case of adiabatic coupling 
the molecular coordinate must coincide with one of 
the normal coordinates q ~ Qy and the respective 
parameter is thus S^y = 1). 

From the general form of system Hamiltonian ([T]) 
one can see that a problem of thermal excitation of 
a molecule in a condensed medium is reduced to the 
problem of description of evolution of harmonic quan- 
tum system ([2]) having the randomly alternated energy 
levels ([3]). By energy conservation law (principle of mi- 
croscopic equilibrium) molecule exhibits persistent ex- 
change of phonons with a heat bath Q through the 
one-phonon mechanism ([5]). Thus, for derivation of the 
master equation it is reasonable to use approach that 
has been recently developed in papers [15, 16] just for 
such types of open quantum systems. 

2. Master equation 

To perform an analysis of behavior of a molecule dur- 
ing the process of its thermal excitation the master 
equation for evolution of observable state occupancies 
7„(t) is required. Since the frequency of transitions 
between the molecular energy levels is the stochastic 
quantity, observable state occupancies appear to be av- 
eraged over random realizations of shifts of the energy 
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levels: 7„(t) = ((r„(i))), where ((...)) denotes the aver- 
aging over random shifts, while r„(t) = {n\p{t)\n) repre- 
sent the non-averaged state occupancies. The molecular 
nonequilibrium density matrix p{t) — trTPs+xit) is de- 
termined as a trace over the thermal bath states (trT) 
on the nonequilibrium density matrix of the entire sys- 
tem ps+xit). In turn, evolution of the entire system is 
governed by the Liouville equation: 



PS+rit) = -iC{t)ps+T{t) 



(7) 



where £(t) = {l/h)[H{t)^ ..] is the Liouville stochastic 
superoperator that acts in the space of both dynamic 
and stochastic variables of the entire system Hamilto- 
nian H{t) ([T]). Expanding ([7]) into diagonal and off- 
diagonal components and solving obtained equations 
for diagonal component of the molecular density ma- 
trix pd{t) = {trTPS+T{t))diag = Z]n (0 1") (" I yields 

Pd{i) = ~ Sldt'trT[fdCvS{t,t')CvfdPs+T{t')], where 
Cv = {l/h)[Vrei,.:] and S{t,t') = e^'' '^^(i-^'')^^ 
with Td being the projection superoperator which trans- 
lates any operator into its diagonal component. 

In many practical cases it is common that a character- 
istic time of transition processes Ttr significantly exceeds 
the characteristic time of thermal relaxation tt- This al- 
lows one to use a factorization TdPs+rit) = pd{t)pT, 
where px = exp{—HT/kT)/trTexp(-'HT/kT) is the 
bath equilibrium density matrix. Thus, in Born approx- 
imation for the Liouville equation one arrives to the fol- 
lowing coarse-grained master equation: 

1 /■* 

Pdi^)=-J^ I dt'trT{Td[VreuU{t,t')[VreUPd{t')pT]x 



xC/+(i,t')]} 



(8) 



where U{t,t') = fexp[-{ilh) J^, dT{Ho{t) + Ht)] and f 
is the Dayson's chronological operator. 

Taking into account (jS]) and using exact form of system 
Hamiltonian and ([5]), one can derive the equation 

for non- averaged state occupancies r„(t): 

f„(t) = / dt'{g {t,t')Tn+l{t')+g n—ln 



exhibit stochastic behavior through the stochastic func- 
tionals J>„„(t,i') — exp{i J^, dT[(m — n)(cj(T) — -cuo)]}, 
while Rx{t) = A^(rjA)e**^^^ + [N{nx) + Ije"**^^^ cor- 
respond to the correlation functions for one-phonon 
transitions between nearest vibrational levels with 
N{n) = [exp{nn/kT) ~ being the Bose distribution 
function. 

The main difficulty in solving of the integro- 
differential equation ([9]) lies in a non-Markovian char- 
acter of the respective integrands. Moreover, one has 
to make explicit averaging of the stochastic function- 
als Tn{t) and J-'nmit,t'). However, when considering 
one-phonon transitions in the second order of pertur- 
bation theory over relaxation interaction, then the non- 
Markoviaty can be neglected [15]. Using this fact and 
taking into account that characteristic time of stochastic 
processes Tgtoch is of the order of thermal relaxation time 
Tt ~ Tstoch and consequently Tstoch «Ttr, one can per- 
form stochastic averaging of © in the non-explicit form, 
that yields: 

^7n(t) + l)T4^+7n+i W + nW^-7«-i W- ^^Q^ 
- [nW+ + {n + l)W4-f^[t) 

The respective averaged probabilities of transitions from 
the first to the zero (+) and from the zero to the first 
(-) vibrational levels are 

W±=2Rey2^xmonx dre^'^"' Rx{t)F{t) (11) 

where F{t) = {{exp{i /J" (iT'[a;(r') — njo]})) is the aver- 
aged correlation function of stochastic shifts. 

As one can see, a general problem reduces to the 
problem of determining the model for levels' stochas- 
tic alternation uj{t), subsequent calculation of the cor- 
relation function F{t) and final evaluation of transition 
probabilities W±. Note that in the case of thermal ex- 
citation, initially only one energy level (zero) is occu- 
pied: 7„(0) = Sno- Thus, evolution of both the mean 
vibrational energy (Ey^bit)) = J2n^'^o1n{t)(n + 1/2) 
and state occupancies 7„(t) are determined by the well 
known expressions [4,5]: 



[Gnn+lit^t') + g„„^lit,t')]rnit')} 



(9) 



Here e„™(t,t') = ^ReJ^xXliin + l)Smn+i + 
n5,„„_i]e'^«(™-")(*-*')i?A(t - f)Tnmit,t') are the 
non-Markovian transition probability densities that 



' {E,^b{t)) = hwo[N{mo){l - e-"*) + 1/2] 



7n{t) 



(12) 



In the above equations k = W+ — W- is the effective 
rate constant of thermal excitation which fully describes 
kinetics of relaxation processes. Therefore, the quantity 
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K represents a major kinetic parameter of interest which 
has been sought as a general solution of the problem 
considered in the one-phonon approximation. It should 
be noted that state occupancies could also be represented 
in the equivalent form 7n(i) = [{ln{t)) — l/2]"[(7ri(i)) + 
l/2]-"-i, where (7„(t)) = N{tuo){1 - e""*) + 1/2 is the 
average mean occupancy vibrational number. 

3. Rate constant 

Let us use symmetrical dichotomous process a[t) as a 
model for stochastic shifts of the molecular frequency 
oj{t) Q from its mean value wq. During such process, 
the frequency exhibits jumps between two equiprobable 
values uj{t) — tuq = a{t) = ±a at random times with 
mean frequency v. It can be shown [15] that in this case 
in (dll) F{t) = (fcie^^r _ j^^^-kiry (^f,^ _ ^^^^ .^^j^^^.^ 

ki,2 = ^(j^± Vi^^ - 40-2). Then, according to ([12]), the 
rate constant can be represented in the following general 
form; 

K :=2tuoCT^!/ VCa^^aIttt^ 5^5- — jtt^ 



(13) 

From one can see, that rate depends on dynamic 
{ujQ, Qxi Ca), shifted {ujq — ujq + a), and stochastic (cr, v) 
system parameters. The last additionally provide for 
rate temperature dependence, as cr = <^{T) and v — 
v{T). Therefore, to perform further analysis, one has to 
represent rate constant in the convenient form. 

First of all, let us make summation over the nor- 
mal modes A in (fO)l . Let the corresponding factor of 
non-adiabaticity, that characterizes the constant of in- 
termolecular anharmonicity ([6]), has a structure 

6 = 2r!A/(r!A) (14) 

where the function /(JIa) = ^/{^\ + "H^) determines a 
density of spectral distribution of frequencies VL\. In 
this distribution the width parameter 77 plays a role of 
adiabatic width for heat bath phonon spectra (in (14) 
the adiabatic limit ^a 1 arrives at rf). Concrete 

value of T] depends on the model chosen for adiabatic 
interaction in given system. If vibrations in molecule and 
medium merge into the unified set of vibrational states, 
then a heat bath with an infinite set of normal modes 
forming an almost continuous spectrum can be chosen as 
a model for the molecular environment. The continuity 



of such spectrum means that after the transformation of 
a sum over A into an integral over and its subsequent 
calculation the value of 77 must be turned to zero: 77 — >■ 
+0. 

This approximation is usually adequate in homoge- 
neous structures of similar molecular nature, for exam- 
ple, in crystals. Another situation could occur in hetero- 
geneous systems (disordered and amorphous structures, 
biological macromolecules etc.), especially at finite tem- 
peratures. In such systems adiabatic spectral widths 
1] — r]x would have additional dependence on the den- 
sity of distribution of normal modes A. In particular, 
their values might coincide with frequencies of certain 
adiabatic (mechanical) modes 77^ — ^^a and would be 
much smaller (77^ << a) than the standard deviation 
for thermal fiuctuations a. On the other hand, values r]x 
should considerably differ from nonadiabatic (optical) 
modes rjx << a << fix- But in all cases, the shifted fre- 
quency of molecular vibrations -cuq = WQ-l-cr, that linearly 
rises with mean random field cr, significantly exceeds the 
adiabatic width: tuq >> rjx- Thus, in the one-phonon 
approximation fix — '^0 the values of 77^ must be con- 
sidered small for all A. 

Note that physically, in the case of uniform distribu- 
tion of normal modes, a width of adiabaticity 7/ forms 
in fact the finest scale over the frequency axis fix- This 
allows one to make in (13), (14) a transformation to the 
variables fix = (Ai7)A = 7;A, where the elementary fre- 
quency shifts AD, = JIa+i — ^\ = V will not depend 
on the number of vibrational mode A. By so doing, 
one can set in (13) the sum over A infinite and, tak- 
ing into account the continuity condition 77 = dfl — > 0, 
make a transformation from a sum of functions fi^x) 
to the integral of a function /(51) by the following rule 
Er=i/(^^A) = liui^^oi^ f^d^]- a resuh, 
with the account of (14) elementary calculations yield: 

4 r^+^ V 

k=-(wo + ct)/ rfz-[(l-z2)2+ (15) 

TT Jo cr cr- 

We see that, since the sum (13) does not generally 
depend on the Bose distribution function for molecular 
vibrations iV(n7o), the very knowledge of particular non- 
adiabaticity parameter S^x (14) allows us to reduce rather 
complicated equation (13) to the simple integral (15), 
which only contains the unknown dimensionless stochas- 
tic factors z = il/cr, vja^ and vj^ja. As a second step, 
let us explore a concrete model for the temperature de- 
pendence of stochastic parameters. For this purpose we 
use thermodynamic approach that was originally pro- 
posed in [15, 16]. In this approach, standard deviation 
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of an amplitude of dichotomous process a is linked to 
the value of mean square of energy fluctuation 6E^ via 
relation: Tr^/i^cr^ — SE'^. On the other hand, from the 
theory of thermodynamic fluctuations in canonical en- 
sembles we know that the latter quantity is given by 
= kT'^{dE/dT), where E = 2TThv is the mean en- 
ergy of dichotomous fluctuations with frequency v. In 
the classical limit, mean energy E — hujc[N{ujc) + 1/2] 
per separate degree of freedom must turn out to thermal 
energy /cT, while in the quantum limit, the same quan- 
tity must correspond to energy of zero-point oscillations 
with some correlation frequency lOc- Therefore, indepen- 
dently of the molecular system analyzed, the stochastic 
fluctuation parameters can be determined by the follow- 
ing expressions 

\<y= [ujJ^)^N{loc)[N{lo,) + 1\ 
\v^{lo,/2tt)[N{uj,) + 1/2] 

General relations and allow one to exam- 
ine important limiting cases. Thus, since a charac- 
teristic scale for stochastic values is determined 
by a frequency lOc- in the case of sufficiently low tem- 
peratures {fujJc » kT), the upper limit of the in- 
tegral can be turned to infinity. In this case 
/„-dzf[(l-zY + ^^']-^=V2,and 

TT 

Conversely, at sufficiently high temperatures {hwc << 
fcr) numerical calculation of the integral in ((T5t yields 

K=l(XA^\)uj^ + Aa)^o(^) (18) 
TT a 

where A ~ 2.10357. Note, that the value of integral in 
(fT5|) rapidly tends to 7r/2 with increasing of TujJq / kT. It is 
equal to 7r/2 up to fourth order even when hujo/kT =^ 1. 
Therefore, practically for all temperatures of interest the 
corresponding thermal excitation rate constant k can be 
well described by quite simple analytic relation (jl7p . 

It is worth noting that equation (|15p can be proved 
only by certain approximations, that simplify general 
character of thermal excitation. Particularly, we use 
both the model of one-phonon transitions between near- 
est energy levels and dichotomous mechanism of fre- 
quency alternation due to the interaction with random 
motions in the environment. 

The first approximation is a direct consequence of bi- 
linearity of the intermolecular anharmonicity operator 
over mutual deviations of the nuclei of molecule under 



consideration and molecules in the environment from 
equilibrium positions ((S]). As it is known, an account 
of non-linear terms significantly complicates the prob- 
lem by making a relaxation of the system essentially 
multi-phonon. In particular, transitions over the one 
vibrational level become possible, as well as transitions 
accompanied by simultaneous creating or annihilating 
of several phonons in the medium [5, 6]. However, tak- 
ing into consideration such effects is important only at 
very high temperatures, when the vibrational structure 
of molecule becomes undefined in a quantum sense. This 
is also the case for molecules with internal anharmonic- 
ity, when only the lowest vibrational levels are well de- 
fined, while energetically higher levels merge into a con- 
tinuous spectrum. 

The second approximation is not critical. Thus, using 
in ([TTjl instead of dichotomous model a common Gaus- 
sian model with analogous in their meaning stochastic 
parameters a and one achieves rate constant in the 
form 

4 /-I+wo/ct 2 

K=-{ujo + a) dz-[l + ^z^]-^ (19) 

TT Jq a 

This equation due to the 6-\ike character of its integrand 
can exactly be reduced to the same form ([T71) as in the 
low temperature limit of dichotomous model ([T5|) . 

In general, the presence of two terms in expressions 
(15), (17), (18) for the rate constant k physically means, 
that in the harmonic (one-phonon) approximation the 
process of molecular thermal excitation is governed by 
the two additive mechanisms. The first is a quan- 
tum mechanism (^ 2a;o) being essentially temperature- 
independent, while the second is quasiclassical one (^ 
2a) being activated with temperature. The mutual con- 
tribution of these mechanisms is determined by charac- 
teristic frequencies of vibrations in the molecule (wq) and 
medium (ujc) with relation to thermal frequency ujt = 
kT/h. Thus, in the quantum limit loq > iOc >> ujt, the 
inequality >> a holds, and the main contribution 
to thermal excitation is introduced by quantum transi- 
tions between nearest vibrational levels of the molecule. 
In this case k = 2Eo/h (17), that in a full accordance 
with the uncertainty principle corresponds to the rate 
speed limit for quantum transitions between orthogonal 
states with energy dispersion Eq — fkoo [17, 18]. In 
the limit of classical transitions ujc < too << a < lot 
and K, = 2ujt/it = 2kT/'!ih, that coincides with the Ein- 
shtein relation for adiabatic rate of the Brownian os- 
cillator relaxation ]13, 14]. In this case, if holding the 
condition lJc < loq, then according to adiabatic theo- 
rem there are no transitions between the energy levels 
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of the molecule. As a result, thermal excitation occurs 
via simple increase of the shifted vibration frequency 
vjQ = Wo + (J « (T = kT/nh due to its direct coupling to 
random motions in the environment. 

We have to note, that quantity k depends on the 
Planck constant h only through the characteristic time 
of thermal relaxation tt = OJ^^ = h/kT involved in 
the Bose distribution function for correlation fluctua- 
tions in the medium N{cOc) = [exp{cOcTT) — This 
property allows one to make a safe transformation from 
the quantum (T — >• 0) to quasiclassical {h — >■ 0) limit 
in a heat bath, considering processes of thermal relax- 
ation to be very fast >> > Wc {tt — > 0) or 
setting temperature relatively high T — >■ oo. By turn, 
the corresponding transformation from the temperature- 
independent to linear in temperature asymptotics for k 
is "well-rcproduccd in the unified manner. Therefore, one 
does not need neither knowledge of the molecular {ojq) 
or correlation (ujc) vibrational frequencies, nor any in- 
troduction of additional factors for "correction" of the 
temperature dependence of rate constant in the classical 
limit [11]. This characteristic feature solves the so called 
harmonic oscillator paradox [19]. Indeed, one may think, 
that in accordance with the quantum theory the overall 
probability of transitions between the nearest levels of 
oscillator, for example between the nth and the (n— l)th, 
should be proportional to the total width of these levels 
(i.e. to the factor ffxlC^n - l)(2|n)(n| + 1)) (5), (11). 
But while turning to the classical limit n — > oo, that is 
achieved at high temperatures (when ojq « kT/h — ^ oo 
or 7V(wo) ^ (|n)(n|) — > oo), this factor diverges. Thus, a 
concept of classical harmonic oscillator as a dynamic sys- 
tem being at the extremely high temperature and hav- 
ing an infinite set of densely distributed levels makes no 
sense. However, from general equations (12) it follows, 
that an evolution of both an average vibrational energy 
and oscillator state occupancies on thermal excitation 
depends on the rate constant n which is determined not 
by the sum, but by the difference between level widths, 
and nor between the some nth and the (n — l)th levels, 
but only between the first and the zero levels at any tem- 
perature and any density of level distribution. Because 
of this, both the numbers of vibrational levels and the 



respective Bose functions for molecular vibrations N{uiq) 
remove from explicit consideration introducing any con- 
tribution to the value of rate constant k (13), that limits 
the process of establishing a final thermal equilibrium in 
the molecule. 
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